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Abstract 


We show that there exists a divergent correlation length in 2d quantum gravity for 
the matter fields close to the critical point provided one uses the invariant geodesic 
distance as the measure of distance. The corresponding reparameterization invariant 
two-point functions satisfy all scaling relations known from the ordinary theory of 
critical phenomena and the KPZ exponents are determined by the power-like fall 
off of these two-point functions. The only difference compared to flat space is the 
appearance of a dynamically generated fractal dimension d^ in the scaling relations. 
We analyze numerically the fractal properties of space-time for Ising and three- 
states Potts model coupled to 2d dimensional quantum gravity using finite size 
scaling as well as small distance scaling of invariant correlation functions. Our 
data are consistent with dh = 4, but we cannot rule out completely the conjecture 
du = — 2«i/a_i, where ct_ n is the gravitational dressing exponent of a spin-less 
primary field of conformal weight (n + 1, n + 1). We compute the moments ( L n ) and 
the loop-length distribution function and show that the fractal properties associated 
with these observables are identical, with good accuracy, to the pure gravity case. 
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1 Introduction 


1.1 General framework 


Two-dimensional gravity has been intensely studied the last six years, both as a 
toy model for four-dimensional gravity and because of its importance in string the¬ 
ory. Let the partition function of two-dimensional quantum gravity coupled to a 
conformal field theory be defined by 

Z{A) = J V[g\D g( f) e - A f d2 ^~ s M(<?>’9\ (i) 

where A denotes the cosmological constant, the integration is over equivalence classes 
of metrics \g\, and Sm( 4>, g ) is the matter Lagrangian. Furthermore, if if is a primary 
conformal field of scaling dimension A 0 in the theory defined by Sm(< f>) in flat space 
and if ((•) ) M denotes the functional average of an observable (•) calculated with the 
action Sm{4>), we have 

M£iM6 )>m ~ _^ 2 |2Ao ~ ( 2 ) 

If V denotes the volume of flat space-time Eq. (|2j) implies the finite size scaling 

j v j v d 2 C ~ r 2 - A » (3) 

for V sufficiently large. In the seminal papers of KPZ and DDK [I], |3|] it was shown 
how to extend this result to the situation where the conformal field theory was 
coupled to two-dimensional quantum gravity. Let Z(V ) be the partition function 
(i, only restricted to universes of space-time volume V, viz. 

Z(V) = Jv[g}V g cj>e- SM ^5(Jd 2 ^-V), (4) 

and let ((■) ) M+G denote the functional average of an observable (•) in the ensemble 
defined by Z(V). With this notation the finite size scaling relation of KPZ and 
DDK reads: 

( f d 2 ZVs f d 2 ('Vs vK)vK')) ~ r 2 ~ A , (5) 

\J J / M+G 

where the exponent A is related to Aq via 


^ a/1 — c + 24 Ao — a/1 — c 

VY^~C 

In (s 6) c denotes the central charge of the conformal field theory defined in flat two- 
dimensional space-time by Furthermore it was shown that the partition 

function for fixed volume behaves as 


Z(V) ~ W 3 , 


7 


c - 1 - yj(c- l)(c- 25) 
12 


(7) 
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We want to emphasize that all these results refer to the partition function or to 
correlators which are integrated over all space-time. In flat space-time the scaling 
properties of the two-point functions (^]) are considered the underlying reason for the 
finite size scaling relations such ([|). Clearly Eqs. (^) and (|7|) are corresponding finite 
size scaling relations in the conformal held theories coupled to quantum gravity and 
it is natural to expect that the scaling properties are dictated by two-point functions 
depending on the geodesic distance. Nevertheless the analogy of @ in quantum 
gravity has only recently been analyzed []|, despite the fact that the two-point 
functions depending on the geodesic distance probe the metric properties of space- 
time in a much more direct way. In fact, as we shall see, the two-point functions are 
perfect probes of the fractal structure of quantum space-time, and they will highlight 
the fact that even if we start out with an underlying two-dimensional manifold, there 
will be metric properties of two-dimensional quantum space-time which cannot be 
viewed as two-dimensional. 

A basic property of the continuum two-point function (Q) (for infinite volume) is 
its invariance under translations and rotations. If the volume V is sufficiently large 
we can write: 




, m-e\-R ) 


n 0 (R; V ) 


where n 0 (R; V ) denotes the volume of a spherical shell of radius R. 
we have from (||) that 


R 2A o ’ ^ 

By definition 


J dR n 0 (R ; V) G^(R- V ) ~ E 1-Ao . (9) 


While it is difficult to generalize Eq. (^) to two-dimensional quantum gravity since 
our physical observables have to be reparameterization invariant, Eq. (|j) has a 
simple translation to quantum gravity: 


G V (R) V ) 



n,+ ft V] 


M+G 


( 10 ) 


In this equation the average is defined via the partition function ([3j) for metrics of 
volume V. D g (£,£') denotes the geodesic distance between £ and £' with respect to 
the metric g, and n g (fi: R ; V) is the volume of a spherical shellf] of geodesic radius 
R with center at £: 


nfit-R-V) = [ d 2 e^5(D g (t,Z')-R). (11) 

We can define the average volume of a spherical shell, or the volume-volume corre¬ 
lator, n(R] V ) by 


n(R-,V) = f d 2 i^n g {i-R-V)\ , (12) 

_ \y j / m+g 

1 Notice that such a spherical shell is not necessarily connected. 
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( 13 ) 


and from dimensional analysis we expect that the analogue of Eq. ([?) will be 

rOO 

/ dR n(R] V) GJR; V) ~ V 1 ~ A , 

Jo 

where A is the KPZ exponent in (||). It is sometimes convenient to use the “unnor¬ 
malized” correlator 

n v (R\ v) = (fj d 2 (VaI (') - «)) ( 14 ) 


rather than G^R] V). In this case we have by definition 

roc 

/ dR n v (R] V) ~ F 1_A , (15) 

Jo 

and for ip = 1 we obtain rii(R ; V) = n(R; V ), justifying the name volume-volume 
correlator for the average value of n 5 (£; R ; V). From the definition of n(R ; V) and 
in accordance with A = 0 for <p — 1, we have 


and if we assume that 



dR m (R-, V ) 


r\j 


V, 


(16) 


n,lR.V) R J - for R<t:V 1/d \ 


(17) 


where we denote dh the intrinsic Hausdorff dimension or the fractal dimension of 
quantum space-time, it is natural from Eq. (|13|) to expect 


a AR-, v )~^fx for , (is) 

for (w) 

Eq. W becomes the quantum gravity version of Eq. 

The main goal of this article is to provide evidence for Eqs. (0)-© and to 
determine d h as precisely as possible. However, at this point let us point to a 
subtlety in the definition of dh- It is possible provide an additional definition of the 
fractal dimension of space-time, starting from ([[]) rather than (f|). In this situation 
we do not keep the space-time volume fixed. Instead we define a “global” fractal 
dimension dn by 

(V) R = (/ d2 (-/g) H ~ R d " (20) 


for R < \ i/d,H _ Xhe average in ( |20l) is over all universes where two marked points 
are separated a geodesic distance R , i.e. calculated from the partition function 


Zr( A) = Jv{g]V g( f>e- A f d 2 ^- s ^Jd 2 ^ J d 2 ?^ 8(D g (£, £') - R). (21) 
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In pictorial terms the calculation of (PU| ) corresponds to the calculation of the average 
length of a random walk which has traveled a distance R. A priori there is no 
reason for d h = d H and it is indeed possible^ to find statistical models where d h ^ 
dff- However, as we shall see, for two-dimensional gravity coupled to conformal 
matter it seems that there is only one fractal dimension, valid at all distances. As 
a consequence, the constant coefficients in Eqs. ([I7|) (p~9|) are independent of the 
volume V. 


1.2 Present status 


Before discussing in more details the model and the methods used, let us briefly 
review what is known and what is conjectured so far. 

Using the so-called transfer matrix formulation [3j it has been possible to calcu¬ 
late Z r(A) in the case of pure gravity. It is remarkably simple J4|: 


Z R { A) = A 3 / 4 


cosh(v // AR) 

sinh 3 (v / Ai?) 


( 22 ) 


From the definition of Z R ( A) it follows that it is related to n\{R ; V) by 

P OO 

Z R {A)= dVe~ AV VZ(V) nRR; V), (23) 

Jo 

and it is easy to prove that 

M R i v ) = R3 fi(x), x = (24) 

where fi(x) can be expressed in terms of generalized hyper-geometric functions and 
/i(0) > 0, fi(x) ~ e _x4/3 for i»l. (25) 


From ( |TT| ) we conclude that dh = 4 in the case of pure gravity, while (p2|) implies 
that dfj = 4. Hence a; is a dimensionless scaling variable. It is useful to write Eq. 
(0) as 

n 1 {R)V)=V 1 ~ 1/dh F 1 (x), (26) 


Fi(x) ~ x dh 1 for x -C 1. (27) 

These equations appear as typical finite size scaling relations, ideally suited for 
numerical simulations. In the case of pure gravity one can calculate additional 
correlators of the so-called gravitational descendents [[23|| . 

The final observable which, in the case of pure gravity, can be calculated from 
first principles is the so-called loop distribution function. If P denotes a point on 
the two-dimensional manifold of volume V, we consider the set of points whose 

2 The so-called multicritical branched polymer models are examples where dh = 2 while dn = rn, 

m — 2 , 3 ,.... 
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geodesic distance from P is equal R. Generically, this set consists of several closed 
loops of various length L t . We take the average over positions P as well as over 
the ensemble of all two-dimensional spherical manifolds of volume V (with weight 
e -s M (<t>,g) if matter fields are included). In this way we get a distribution py{R 1 L ) 
such that p v {R, L)dL measures the average number of loops with length between L 
and L + dL per manifold with volume V. We can now define the moments of L: 

poo 

(L n (R)) v = dL L n p v (R, L). (28) 

Jo 

Notice that 

n 1 (R,V) = (L(R)) v . (29) 

In the case of pure gravity it is possible, using again the transfer matrix, to calculate 
pv(R,L ) in the limit where V —> oo 

Poo (R, L) = p(L/R 2 ), (30) 

where 

p(y) = const, x (p” 5/2 + ^ y~ x/ 2 + yp 1/2 ) z~ v ■ (31) 

One observes the singularity of p for small y. This is why a cut-off e in L is needed. 
In fact, (L°) and (L) are singular. If we cut off the integral in (^) at e (the lattice 
length unit of L) we obtain: 


(L(R))v =oo = —i? 3 + 0(l), (32) 

(L n (R)) v=0 o = c n R n , n > 1.. (33) 

If the volume V is finite ( |33| ) is replaced by 

(L n (R))y = V 2n ^F n (x), (34) 

where dh = 4 and where 

F n {x) ~ x 2n for i«l. (35) 

Since we know that L has dimension of length and dh = 4, i.e. R has dimension [ L 1 / 2 ], 
relation (^) is dimensionally correct, but it is clear that we need a dimensionful 
parameter in the relation ( L(R )) ~ R dh ~ l . It can only be provided by the cut-off e. 

It is natural to conjecture |5],|6j that finite size scaling relations like (^6|) and ( |27|) 
are valid also for unitary conformal field theories coupled to quantum gravity, except 
that dh could be a function of the central charge c of the conformal field theory. 
Furthermore, it has been conjectured that the same dh(c) and the same scaling 
variable x = R./V ] ! dh ^ govern the finite size scaling relations of all the correlators 
Gtp(R ; V) of primary conformal fields <p of the given conformal field theory: 

G V (R-, V) = J^KSM = r- a 9„W, (36) 
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where the scaling function g v (x) then behaves as 

g<p(x) ~ x~ dhA for i«l, (37) 

and 

g v {x) —■> 0 for x —> cx). (38) 

It is sometimes convenient to use the “unnormalized” correlation functions instead: 


n lf (R- 1 V) = 


^±F v (x) = V'-^F v {x) 


(39) 


where 


F v (x) ~ x dh 1 dhA for i«l. (40) 

A priori it is not clear what to expect for pv(R , L) if conformal matter is coupled 
to quantum gravity. Rather surprisingly, it seems that the finite size scaling relation 
(|35|) is still valid with dh — 4 replaced by dh(c), and that p(y) is still a function only 
of L/R 2 when V —> oo. In case dh(c) is different from 4 this implies that L has an 
anomalous scaling dimension relative to V. This has been shown convincingly to be 
true for the c = —2 model in ||. In this paper we show that the numerical data 
is consistent with Eq. ([h|) and Eq. ( |35| ) by extracting consistent values for dh from 
(|34l) and by showing that (^) holds with very good accuracy. We should warn the 
reader though, that since dh ~ 4 in our case, the scaling fl34|)-(|35D is an assumption 
based on the c = —2 result. 

It should be emphasized that the assumption of a single new parameter c4 (c) 
which determines uniquely all scaling in the given conformal held theory coupled to 
quantum gravity is a strong assumption which so far has avoided a rigorous mathe¬ 
matical proof. There has very few attempts to discuss analytically the appearance of 
a divergent correlation length when, say, a spin system coupled to gravity becomes 
critical 0- However, the scaling has already been tested numerically for the Ising 
model and the three-states Potts model coupled to gravity mm- Recall that 
the unitary conformal held theories with central charge c between 0 and 1 are the 
so-called minimal (q,p) rational conformal held theories with (q,p) = (m,m + 1), 
m = 2,3,... and the corresponding central charge is c = 1 — 6 /m(m + 1). Pure 
gravity corresponds to m = 2, i.e. c = 0, the Ising model corresponds to m = 2, i.e. 
c = 1/2 and the three-states Potts model corresponds to m = 5 and c = 4/5. Using 
ni(R; V ) and n^R] V ) it was found that dh ~ 4 independent of these three values of 
c, while dhA was found to be consistent with the value predicted by KPZ, provided 
dh was attributed a value close to 4 (TTJ . In this way the numerical simulations have 
so far provided some support for the scaling hypothesis outlined above, but with the 
rather surprising result that dh ~ 4 for 0 < c < 1. This is in contradiction with the 
two theoretical predictions for dh which exist so far: 


and 


j(n) 

a h 


df = 2 x 
h a/25 


24 


V25^c +V49^c 


C + a/1 — c 


(vT=c + 


(=-)■ 


(41) 

(42) 
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The first formula |TB|, |TT| is derived under certain assumptions about the scaling of 


the diffusion equation in Liouvillc theory and predicts a rather slow change in dh as 
a function of c. The second formula, derived for (q,p) models coupled to quantum 
gravity by means of string held theory |I8|], is based on the assumption that the 
“time” in string held theory can be identified with geodesic distance, ft is true by 
construction for c = 0, but might not be true for c ^ 0. Note in particular that 


d^\c) —> oo for c —> 1 and d^ l \c) —> 0 for c —> — oo. In the table below we have 
summarized the predictions made so far, including numerical results too. 


j(“) 


c 

d 

a h 

d (n) 

a h 

numerical 

— oo 

2 

0 

- 

-2 

3.562 

2 

3.58 ±0.04 

0 

4 

4 

~ 4 

1/2 

4.21 

6 

~ 4 

4/5 

4.42 

10 

~ 4 

1 

4.52 

oo 

~ 4 


Let us emphasize the numerical results for c = —2 [D| (see also |13| for earlier mea¬ 
surements). For c = —2 one can avoid the use of Monte Carlo simulations since 
there exists a recursive algorithm which directly generates independent configura¬ 
tions. The statistics are for this reason very good and the numerical results clearly 
favor the prediction d^ . In addition, the measurement of pv(R,L) for c = — 2 sup¬ 
ports (|34D and (|35]) with c4,(c) = d^\c) and L/R 2 as dimensionless scaling variable. 


1.3 Outline 

It is the aim of the present article to discuss the current available precision in the 
determination of dh(c ) by the use of numerical simulations and to test the scaling 
hypothesis ([H)-([10)- We perform simulations for the Ising model coupled to gravity 
and the three-states Potts model coupled to gravity as well as simulations for pure 
gravity (as a check of the accuracy). In particular we measure the volume-volume 
correlator ni(R]V), the moments ( L n (r))v (n = 2,3,4) as well as the spin-spin 
correlators G<p{R] V) and n^R; V) and apply finite size and small distance scaling in 
order to determine c4,(c). The existence of such scaling provides evidence that there 
is a unique fractal dimension d^(c) which extends over all distance, a highly non¬ 
trivial fact for the critical spin systems coupled to quantum gravity. Furthermore 
we measure the so-called loop length distribution function pv(R,L) which provides 
additional insight in the fractal structure of quantum gravity. By comparing the 
different measurements (which are more or less consistent) we get a handle on the 
systematic finite size effects which affect the determination of dh(c). 

In the next section we introduce our discretized model of quantum gravity, the so- 
called dynamical triangulation model of quantum gravity, and discuss how to extract 
the observables defined above within the framework of dynamical triangulations. 
In particular we discuss the so-called shift a which appears when the observables 



















measured on a finite lattice is mapped to the continuum observables defined above. 
In section 3 we describe the numerical method and present our results. Finally we 
discuss our results in section 4. 


2 The model 


We define the regularized theory of two-dimensional quantum gravity via the for¬ 
malism known as dynamical triangulations ||. Some of the exact results for pure 
gravity mentioned in the introduction were derived within this formalism which in 
certain respect is more powerful than a continuum formalism, even for analytical 
calculations. 

In this formalism surfaces are constructed from equilateral triangles glued to 
together to form triangulations with spherical topology^. We allow for the generation 
of certain degenerate triangulations. They are described most easily by using the 
fact that the graphs dual to triangulations are 0 3 -graphs. Regular triangulations 
correspond to connected </> 3 -graphs without tadpole and self energy sub-graphs. The 
degenerate triangulations include the connected 0 3 graphs with tadpoles and self 
energies. The use of the full set of connected </> 3 -graphs of spherical topology is 
known to improve scaling compared to the class of 0 3 -graphs which corresponds 
to regular triangulations |22j]. The coupling of spin systems to gravity is done by 
assigning a spin to each vertex in the triangulation, either Ising of three-states 
Potts spin, depending on the model. At the critical temperatures for the spin 
systems on dynamical triangulations such models describe conformal held theories 
of central charge c = 1/2 and c = 4/5 coupled to gravity, respectively. In this way 
the discretized partition function corresponding to the formal continuum partition 
function (|) is 

Z(N) e~ ST » ({<7 ' } \ (43) 

TN {cTi} 

where the first summation is over all (spherical) triangulations Tjy consisting of N 
triangles and where the second summation is over all spin configurations on T/y. 
At ({A/}) denotes the action for the Ising or the three-states Potts model defined on 
the triangulation T in the way mentioned above, viz. 


<St({w}) = P a k°u 

(k,l) 


(44) 


where f3 is the temperature, the summation is over all pairs of neighboring vertices 
(k,l) and q refers to the g-states Potts model (q — 2 for the Ising model). For the 
g-states Potts model the spin variables can take q values and oyoy is a symbolic 
notation for 

(45) 


OfcCp = 


Q 


1 


q-1 ak,ai q-1 
The continuum volume V is related the number of triangles N by 

V oc Ne 2 


3 We will consider only surfaces with spherical topology in the following. 
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where £ is the length of the individual links in the triangulations. As usual in lattice 
theories we work in units of e, i.e. £ — 1. Phased entirely in the framework of 
statistical mechanics we study the annealed average of g-states Potts models on the 
class of all triangulations with spherical topology. 

Geodesic distances r on the triangulations are defined either as the shortest link 
distance between two vertices or as the shortest path through neighboring triangles. 
While these two distances can vary a lot for two specifically chosen vertices in a 
given triangulation, they are proportional when the average is performed over the 
ensemble of triangulations. We will report here the results obtained by the use of 
link distance, since it is known that the short distance behavior of correlators suffer 
from less discretization effects if we use the link distance than if we use the triangle 
distance^ || . 

The lattice analogue of the correlator functions ni(R ; V), n v (R] V ) and g^{R] V ) 
can now be defined. We will use the same notation for the function names, but r and 
N instead of the continuum variables R and V. No confusion should be possible. 
Let i and j denote lattice vertices and the (lattice) geodesic distance between i 
and j and r an integer (the lattice analogy of R, measured in units of e). Then, 


7ii(r; N) 
n<p(r; N ) 
G^t-N) 



(46) 

(47) 

(48) 


where cycr, is defined as in ([15]). For a given vertex i we always perform the sum¬ 
mation J2j in Eqs. (|46|) - ([f8|) . ffowever, from a numerical point of view it is not 
convenient to perform the summation over i for a given triangulation since the sums 
over j for different choices of i will be highly dependent. On the other hand it is 
not convenient either from the point of view of numerical efficiency to choose just a 
single vertex i for each triangulation. We found it convenient to generate indepen¬ 
dent configurations (by Monte Carlo simulations) and for each such configuration 
to choose randomly a suitable number Nq of vertices i and replace 


1 

N 


N 


E 


N 0 


— y 

N 0 ^ 


a=l 


We have checked that the numerical results are idependent of No . In particular, one 
can choose the extreme value N 0 — 1, although it is not convenient from the point 
of view of efficiency, as mentioned above. 

4 Intuitively it is to be expected that the link distance will behave more “continuum like” for 
short link distances compared to short distance triangle distances. At short distance the triangle 
distance is quite “rigid”: each triangle has 3 neighboring triangles (except for some degenerate 
triangulations), while vertices can have a variable number of neighboring vertices. 
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We now expect the following short distance behavior: 


TT-i (r; N) ~ r dh 1 , 


n v {r\ N) ~ 


G v (r;iV)~r-^ A , 


(49) 


as long as r -C N 1 ^^, where the constant coefficients in the above relations are 
independent of N. The scaling hypothesis can be formulated as 


where 


ni(r; N) 

~ N'-V^F^x), 

(50) 

n v (r- N ) 

~ N'-^FF^x), 

(51) 

G^r-N) 

~ N~ A g v (x), 

(52) 


x - JV*/*.' 

(53) 


and the functions F(x), F lf (x) and g v {x) should approach their continuum counter 
parts in the limit N —> cx). ft is known that this approach can be slow |5|.[6[]. It 
can be improved by applying a typical finite size scaling argument 0 : we can only 
expect strict proportionality 


re 2 l dh oc R and Ne 2 oc V 


(54) 


for N —» oo. For finite N one expects 

R 

= const, x 


yVd h 


N^/dh 




(55) 


since 1/N 1 Fh measures the linear extension of the system, i.e. the simplest finite 
size correction to x can be parameterized as 


r + a 

x = ’ 

where a is called the “shift”. In principle this shift might depend on the observable 
we consider. This shift was first introduced in [ 6 | in the case where r is the triangle 
distance. It was shown that in the case of pure gravity a remarkable agreement 
between the discrete and continuum volume-volume correlators can be achieved for 
a ~ 5.5, even for quite small lattices. These authors introduce a second shift variable 
N 1 F h —> N 1 / d h _|_ 5 which in this study we will set equal to 0. In the case of link 
distance, as our results will show, a is much smaller, in agreement with the remark 
above concerning finite size effects for link and triangle distances, and it is possible 
to obtain very good results from finite size scaling without introducing it |5|, pT|p!2[| . 
Its role, however, is quite important when trying to determine d^. One cannot 
ignore it since, as it was shown in |J, it is a necessary finite size correction in the 
case where one uses triangles to define geodesic distance and area in the discretized 
version of ni(R] V). Its introduction improves the fits in the case of link distance 
as well, but it will add an extra free parameter, since we have no analytical results 
in the case of matter coupled to gravity. This will make the determination of dh 
more difficult. Moreover, a improves dramatically the small distance scaling of the 
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correlators, making it possible to perform direct fits near the origin in order to 
extract the scaling exponents, which would otherwise be impossible. 

Finally, the continuum definitions for the loop length distribution, Eqs. (30)-(33), 
are valid at the discretized level with the replacement 


R —> const, x r, L —> const, x I, e = 1, 


(57) 


where r is the lattice geodesic distance (the link distance) and l is the lattice length 
version of L. It is defined in the following way: for a given triangulation let j be a 
given vertex and let ji, ■ ■ ■ ,j n be the set of vertices located a link distance r from 
vertex i. The vertices j i,..., j n can be divided in connected maximal subsets such 
that j a belongs to the subset {jp 1 ,..., jp k } if it is neighbor to any of the jp s and 
not neighbor to a j a in {j\ ,..., j n }/{jg l , • • •, jp k }■ For such a maximal subset we 
define l = k. The loop length distribution pAr(r,/) at the discretized level measures 
the average number of such connected boundaries of length / corresponding to the 
geodesic radius r for the ensemble of triangulations with N triangles. Corresponding 
to (^8D we have 


<!»>» = E^W), 

l 

(58) 

while the scaling hypothesis (^) becomes 


(r (r)) N — N 2 n / dh F n (x), 

(59) 

where 


F n (x) ~ x 2n for i«l. 

(60) 


As already mentioned, we will provide evidence that the loop length distribution in 
the limit N —> oo, even after gravity is coupled to matter, has the form 


pN^oir, o = 4 /%)> y = 4- ( 61 ) 


Numerical simulations of pure gravity give good agreement with p(y) given by m 
Here we will be interested in a determination of py(R,L) for the critical Ising 


19 


and the critical three-states Potts model coupled to quantum gravity. 


3 Numerical Simulation 


The numerical simulations are performed as follows: The Monte Carlo updating 
of the triangulations is performed by the so-called flip algorithm and the spins 
are updated by standard cluster algorithms. The flips are organised in “sweeps” 
which consist of approximately Nl accepted flips where Nl is the number of links 
of the triangulated surface. After a sweep we update the spin system. All this is 
by now standard and we refer to |2l],[22j for details about the actions or Monte 
Carlo procedures. We use the high quality random number generator RANLUX |29| 


whose excellent statistical properties are due to its close relation to the Kolmogorov 
K-system originally proposed by Savvidy et.al. 


28 in 1986. 
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All our runs for the Ising and three-states Potts model were made at the exactly 

urn 


known infinite volume critical temperatures /3 C 


The sizes and the number 


of sweeps we use are different depending on the observables that we measure. The 
largest amount of statistics we gathered were for measuring ni(r\N) and n^(r; N) 
for the Ising and three-states Potts models coupled to gravity, were we performed 
1.7 — 5.0 x 10 6 sweeps on surfaces with 16000-128000 triangles. We also report 
results obtained on surfaces with 256000 triangles with a smaller number of sweeps 
(~ 0.8 x 10 6 ). For the moments ( l n ) we needed much less statistics: We performed 
3.0 — 6.0 x 10 5 sweeps for the 16000-64000 lattices and 3.5 — 7.0 x 10 4 for the 128000 
lattices. Although less statistics was necessary in the latter case, the computer effort 
needed for the measurements is quite significant compared to that of measuring 
ni(r; N) and n ¥ ,(r; N), especially for the largest lattices. 

We can extract the Hausdorff dimension from the scaling hypothesis 


n,(r;JV) = n v (r;N)=N 


1—A— 


,/dk F v (x) 


(62) 


as was done in [pj.jB|]. The analysis of Eq. (^) is performed by “collapsing” n-\_ tV (r; N ) 
for a given number of lattice sizes. A fit to p n (x)e~ ax , thought of as an interpolating 
function, is performed for given ( dh,a ) and a value for the x 2 (dh,a ) is obtained. 
p n (x) is a polynomial of x of order n. n is chosen large enough to capture the 
functional form of N) by checking that X 2 (c4,a) does not depend on n for 

a range of n and small enough in order to leave enough degrees of freedom. The 
errors which we use in the determination of y 2 (d/ l ,a) are computed by binning our 
data. We refer the reader to the figures in ||,^] in order to appreciate pictorially 
the impressively good scaling that ni iV ,(r; N) exhibit in the simulations. 

In Fig. [I] and Fig. [2| we show the results from the analysis of ni )V ,(r; N). We 
see that determining a is crucial for extracting d h . In these figures we show c4(a) 
given by the value of dh. which minimises X 2 {dh,a) for fixed a. The errors are 
computed by the interval of dh which changes x 2 ( a ) max{2,2y; 2 (a)} where 

X 2 ( a ) — min d h {x 2 (dh,a)}- As we can see, the value of dh changes considerably 
with a so we need to determine the range of its acceptable values, since a = 0 is 
by no means a special choice in our consideration. We do this by minimising % 2 (a) 
which gives the results shown in Table |T| and Table ^]. The errors quoted are com¬ 
puted by considering the interval of a which changes Xmin max{2, 2 Xmin} where 
Xmin = m in a {x 2 (a)}. Notice that a = 0 is not significantly far from the optimal 
choice of a. This explains the good quality of the results reported in [^,|l]J in the 
case of link distance. This is not true if one uses the triangle distance for the lattice 
sizes we consider here (we expect a to become numerically less significant for large 
linear size of the system). 

dh can also be determined from the small x behaviour 


F^x) 


x 


dh~ 1 


F v (x) 


x 


d h ( 1-A)-1 


(63) 


and 


d log F x (x) 
d log a; 


dh — 1 + x k + ... , 


(64a) 
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dlogF^x) 
d logx 


— dh( 1 — A) — 1 + x l + 


(64b) 


The calculation of the logarithmic derivative is performed by using a 9 point Savitzky- 
Golay smoothing filter |30| with a 8 th order interpolating polynomial. The errors 
are computed by binning our data. The use of the filter improves the computation 
of the derivatives, especially near the origin. The analysis of the logarithmic deriva¬ 
tives provides an excellent pictorial way for realizing the scaling given by Eq. (|63|) 
near the origin. Fig. shows no such scaling if a = 0.0. Similar lack of scaling 
for a = 0.0 is observed for all correlation functions we analyzed near the origin. 
There is an optimal value of a, however, where scaling becomes manifest. This is 
shown in Fig. |3| and Fig. [4]. In order to determine the optimal a we use direct fits 
to Eq. (^) . The results are shown in Table |5] and Table m In order for the corre¬ 
lation functions to have a “smooth” continuum limit j27| it is very important that 
the value of dh extracted from those fits is the same as the one extracted from the 
analysis of Eq. (|62|) since in principle the two values can be different. For ni(r; N ), 
dh extracted from finite size scaling variable x corresponds to dn of Eq. ([20j) and dh 
extracted from the small distance behaviour corresponds to dh of Eq. (|T7|). In the 
case of ?r ¥ ,(r; N), agreement between the two values implies also the existence of a 
diverging correlation length for the matter fields OH. a fundamental assumption 
for scaling in critical phenomena. For matter systems coupled to gravity we could 
find the situation where, despite the fact that we have a continuous phase transi¬ 
tion, the scale associated with the geometry diverges whereas the scale associated 
with matter does not, as is the case of many Ising spins coupled to two-dimensional 
quantum gravity [I4|]. From our results we conclude that we do obtain the same 
scaling exponents in Eq. (|62|) and Eq. (|63D, the small differences being attributed to 
finite size effects. We also observe that by introducing the shift a, scaling at small x 
with reasonable values for dh appears from the analysis of much smaller lattices than 
it was thought it would be necessary before. Our results on the largest lattices are, 
however, necessary in order to gain confidence that the fits are stable with respect to 
changing the points that one includes in the fits, dh decreases slightly by removing 
points from the origin, but the fits are more stable for the largest lattice that we 
use. 

At this point we would like to present our results for the ordinary spin-spin 
correlation function G v (r]N). It is a pleasent surprise that it exhibits excellent 
scaling properties, in fact better than the “unnormalized” function n^ir^N). The 
reason for the improved quality is presumable that some correlated fluctuations in 
spin and geometry are cancelled in G^r; N ). Our results for the value of d h obtained 
from n v {r\N) are confirmed and give us further confidence on the existence of a 
diverging correlation length for the matter fields. Recall that we expect the scaling 


G r (r-,N)~N~ A g v (x), 


(65) 


where 


w( 


x 


x 


-dh A 


I«l. 


( 66 ) 


5 We should warn the reader that the values of y 2 reported for all fits are to be strictly compared 
only for same lattice sizes since the amount of statistics varies between lattice sizes. 
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Iii Fig. [| we see that the scaling (|(75]) holds very well if one chooses the appropriate 
value of a. g v (x) has a stronger dependence on a near the origin than Fi >ip (x) and 
we do not obtain good scaling at a = 0. The value of a chosen in the plots is taken 
from the fits to Eq. (j66|) . The results of the latter are shown in Table [7]. As one can 
see from the plot of the logarithmic derivative of g v (x) in Fig. (|j), the fits are stable 
over a wider range than those of Eq. (j63|) for n v (r-,N) and we get and excellent 
agreement with the value of dh we have obtained so far. In Table |7| we show our 
results for different cuts for the range in r: We compare the same range used in 
Table |i| where applicable, a typical point in the region of stability and finally the 
point where y 2 becomes of order 1. This result further supports the statement that 
the correlation length for the spin-spin correlation function diverges at the critical 
temperature. 

The results for the Hausdorff dimension are further tested by measuring the loop 
length distribution p N (r,l) and calculating 


</»> 


where 


n ~ N d h F n (x ), n = 2,3,4, 

(67) 

F n (x) ~ x 2n , i«l. 

(68) 


We have analyzed Eq. (|67]) in a similar way we have analyzed ni iV5 (r;IV) by col¬ 
lapsing the distributions. In Fig. [7] we show the collapsed ( l 2 {r)) N function and we 
observe that scaling holds very convincingly. Quite similar plots can be obtained for 
(/ 3,4 (r)) 7 v- We observe that the extracted dh(a ) has a very weak dependence on the 
shift a contrary to what we found for (Z 1 (r))jv before. This is shown in Fig. § for the 
Ising model. Similar graphs can be obtained for all models and moments and we 
show representative values of (4(a) in Table |3|. In Table f| we show the best value 
for the shift a and the corresponding value c4(a). This is slightly different than the 
procedure we followed in the analysis of ni^{r\ N), since now d h (a) is not invertible. 
In Fig. [7] we show that the scaling of Eq. |SD| holds very well near the r origin. 

Finally in Fig. - Fig. [L2| we show our results for the loop-length distribu¬ 
tion function pn(x,1). In Fig. 0 we show the pure gravity measurements together 
with a fit to Eq. ( |3T1) . In the fit we simply rescale y and pN(r,l)r 2 , but we oth¬ 
erwise keep the same coefficients as in Eq. ([H]). The data is consistent with the 
theoretical prediction, as was found before in |L9|, but unfortunately it has no pre¬ 
dictive power to determine convincingly the terms in Eq. (|3ll) for c ^ 0. We are 
able to check, however, some of the scaling properties of pjv(r, Z)r 2 : The part of 
the curve corresponding to the continuum behaviour is independent of N and that 
Pn{t,1){t/N l / dh ) 2 = N°F p (y, L/N 2 / dh ) for all values of y and L/N 2 / dh . This is 
shown in Fig. for the three-states Potts model. It is also curious that in the 
range in y where we observe continuum behaviour, pAr(r, l)r 2 seems to be indepen¬ 


dent of m — 2, 3, 5. In this range the data points in Fig. |lT] are within statistical 
error on top of each other. One could be tempted to conjecture that the fractal 
properties of space-time for unitary theories with 0 < c < 1 are independent of c. 
In the same figure, however, we observe that the finite size corrections do depent on 
c and it is not so clear which is the borderline in the range in y where we observe 
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continuum physics and where finite size effects are important. In contrast, for the 
c = —2 model simulated in [Q[ 10 ] we see a clear difference from our data, reflecting 
the fact that in this model the fractal dimension of space-time is 3.58(4). In Fig. [l^ 
we compare our data with that of P,|T0 . 


4 Discussion 


In this work we have measured the Hausdorff dimension dh for 0 < c < 1 matter 
coupled to two-dimensional quantum gravity using various scaling arguments. In 
the spirit of | 6 ), we introduced the shift a in the investigation, which as a finite 
size correction, improves dramatically the scaling of correlation functions for small 
geodesic distances r giving consistent results for dh■ For finite size scaling, the 
shift must be included in the analysis, even when we use the link distance in the 
definition of correlation functions, and can be used to estimate the systematic errors 
introduced by finite size effects. Its effect is greatest in the case of the two point 
functions N) and G v (r] N ). By studying the scaling behaviour of the spin- 

spin correlation functions we verified and extended the results of [[□], [L2| on the 
existence of a diverging correlation length for matter. We also verified the scaling 
properties of the moments (L n ) for pure gravity predicted in ||,[[| and found that 
similar scaling holds for 0 < c < 1 matter coupled to gravity. We measured the 
loop-length distribution function p^ir,!) and showed that in the continuum limit 
Poo(i?, L) is not or very little affected by the back reaction of 0 < c < 1 matter to 
gravity. 

The results of our measurements of dh are consistent with the earlier observations 
in that the presence of 0 < c < 1 matter has no or very small effect on d^. 
This is in contradiction with the analytic result dh = 2m, as was mentioned in the 
introduction. It has been argued that the reason for the apparent contradiction 
is that a large Hausdorff dimension dh = 2m, m = 3,5,... implies a very small 
linear extension N 1 ! dh for the range of N accessible in the numerical simulations. 
However, if this argument was correct it would be very hard to understand how 
one can measure with excellent precision the correct KPZ exponents for the Ising 
and three-states Potts models from the integrated correlators. Moreover the correct 
scaling of the correlation functions defined in terms of geodesic distance provide 
even stronger evidence that we see the correct coupling of matter to the geometry 
in the simulations, including the fractal properties of the metric. If the linear size 
of the systems were much too small one should not be able to observe the correct 
critical behavior of conformal matter fields coupled to quantum gravity. Moreover, 
as it was already mentioned, it is known [9j.|l0],[l3j that dh = 2m is inconsistent with 
numerical simulations on the c = —2 model coupled to gravity, which corresponds 
to m — 1, i.e. dh = 2. In this case the simulated systems have a quite large linear 
size. 

As it was mentioned in the introduction the alternative prediction coming from 
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scaling arguments for the diffusion equation in Liouville theory, 


d h = 2 x 


v/25^c + V49 - c 
v /25^+ x /T^ ’ 


(69) 


is in excellent agreement with the c = —2 simulations, and it agrees with the rigor¬ 
ously established result dh = 4 for c = 0. 

Our simulations are consistent with the dh = 4 conjecture hif /or Ising and 
three-states Potts model, especially when looking at correlators not involving the 
matter fields. Although the reader can draw her/his own conclusions from our 
measurements, we included a summary of our results in Table | 8 ] where we display 
the most probable range for dh as measured by the different scaling arguments we 
used in this article. We observe that the values of dh coming from correlation 
functions involving matter fields are consistently higher. The errors are big enough 
to include dh = 4 within 1-1.5cr and one can argue that finite size effects are bigger 
in this case: As pointed out in |5[, the scaling behaviour of the spin-spin correlation 
function N ) is the difference fi(r) — 1 /(q — l)/ 2 (r) of the correlators f\{r) and 

f 2 (r ) which count the number of like and different spins at distance r. They both 
scale identically as n\{r\N) and the scaling behaviour of n^{r]N) is obtained by 
exact cancellation of the leading terms. In favour with the dh = 4 conjecture is our 
result that p 7 v(r,/) seems to be independent of m in the region where it exhibits 
continuum behaviour. 

One, however, cannot claim that our simulations exclude the prediction of Eq. &■ 
Except for nfiv, N), one can see that dh increases slightly with m, although the sig¬ 
nal is not clear enough to allow for a clear distinction. There exist the possibility 
that the larger values of dh obtained from n^r; N ) and G^(r; N) are a true signal. 
Firstly, it is of course interesting that the central values of dh actually agree with Eq. 
(|69|), although with large error bars. Secondly, one cannot exclude with certainty 
the possibility that simulations on larger lattices will not shift slightly the values of 
dh extracted from the other observables to the ones of Eq. (|69|). In this connection 
one should bear in mind that the finite size effects for the Ising and the three-states 
Potts model are larger than for for pure gravity, as shown in [|24|. Nevertheless one 
can hope that future simulations might decrease systematic errors to the point that 
one will be able to reject convincingly at least one of the above predictions. 

As we noted in the introduction, the Hausdorff dimension comes from the cutoff 
dependence as L — > 0 of the integral (|32|). The cutoff dependence disappears in 
the integrals (|33|) for n > 1 and one obtains the scaling ([34])-([15]) which seems very 
well satisfied also in the case of the Ising and three-states Potts model coupled to 
gravity, as well as for c = —2! Therefore, in order to be able to fully understand the 
concept of the Hausdorff dimension for c / 0, one needs to provide an explanation of 
expressions like Eq.(^C]) and (pl|) for the loop length distribution function p v (R,L), 
but with different powers of y, also for c / 0. 
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m 

dh 

a 

Nt 

2 

4.05(8) 

0.60(20) 

8000 

64000 

3 

4.11(10) 

0.48(28) 

16000 

128000 

5 

4.01(9) 

0.15(26) 

16000 

128000 

2 

4.05(15) 

0.70(60) 

16000 

64000 

3 

4.11(11) 

0.50(40) 

32000 

128000 

5 

3.98(15) 

0.00(60) 

32000 

128000 


Table 1: The Hausdorff dimension dh as determined from collapsing the n\{r\N) 
correlation functions for pure gravity (m = 2), Ising (m = 3) and three-states Potts 
model (m = 5) coupled to gravity. 


m 

dh 

a 

N t 

3 

4.28(17) 

0.60(30) 

8000 

128000 

5 

4.46(33) 

0.53(51) 

8000 

128000 

3 

4.26(26) 

0.56(48) 

16000 

128000 

5 

4.45(40) 

0.55(95) 

16000 

128000 


Table 2: The Hausdorff dimension dh as determined from collapsing the n^(r;iV) 
correlation functions for pure gravity Ising (m = 3) and three-states Potts model 
(m — 5) coupled to gravity. 



dh 

n 

m 

= 2 

m = 

3 

m 

= 5 


a = 0.00 

a = 0.25 

a = 0.00 

a = 0.25 

a = 0.00 

a = 0.25 

2 

3.88(3) 

3.90(3) 

3.99(4) 

4.01(4) 

4.07(3) 

4.10(3) 

3 

3.94(3) 

3.96(3) 

4.06(4) 

4.09(4) 

4.15(4) 

4.17(4) 

4 

3.95(3) 

3.97(3) 

4.08(4) 

4.10(4) 

4.16(4) 

4.18(4) 


Table 3: The Hausdorff dimension c4(a) as determined from collapsing the ( l n (r)) N 
distributions for N T = 16000, 32000 and 64000 for pure gravity (m = 2), Ising 
(m = 3) and three-states Potts model (m = 5) coupled to gravity. 
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n 

m 

= 2 

m - 

= 3 

m 

= 5 


dh 

a 

dh 

a 

dh 

a 

N t = 16000-64000 

2 

3.88(3) 

0.00(15) 

4.00(4) 

0.0(3) 

4.08(3) 

0.0(3) 

3 

3.94(3) 

0 .0(2) 

4.08(4) 

0.0(4) 

4.16(4) 

0.1(3) 

4 

3.95(3) 

0.00(15) 

4.10(4) 

0.0(5) 

4.17(4) 

0.1(4) 



hP- 

O 

O 

o 

-64000 



2 

3.86(2) 

0 .00(10) 

3.98(4) 

0.05(10) 

4.07(4) 

0 .12(10) 

3 

3.92(2) 

0.05(10) 

4.05(4) 

0.07(18) 

4.16(3) 

0.15(15) 

4 

3.92(2) 

0.05(10) 

4.06(5) 

0.10(18) 

4.17(3) 

0.13(15) 


Table 4: The Hausdorff dimension dh(a) as 
correlation functions for pure gravity (m 
Potts model (m = 5) coupled to gravity, 
corresponding dh(a) is recorded. 


determined from collapsing the (P(r))jv 
= 2), Ising (m = 3) and three-states 
The best value for the shift a and the 
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m 

N t 

dh 

a 

C 

x 2 

T'min 

T max 

3 

256000 

4.081(4) 

0.514(4) 

1.049(9) 

3.7 

1 

6 


128000 

4.098(7) 

0.529(5) 

1 . 01 ( 1 ) 

2.1 

1 

5 


64000 

4.080(5) 

0.517(4) 

1.06(1) 

6.1 

1 

4 


32000 

4.041(6) 

0.492(5) 

1 . 11 ( 1 ) 

7.4 

1 

4 


16000 

3.969(9) 

0.448(7) 

1.25(2) 

8.6 

1 

4 


256000 

4.028(7) 

0.447(9) 

1.19(2) 

1.3 

2 

7 


128000 

3.994(7) 

0.417(9) 

1.27(2) 

1.8 

2 

6 


64000 

3.961(8) 

0.39(1) 

1.36(3) 

1.7 

2 

5 


32000 

3.86(1) 

0.30(1) 

1.65(4) 

5.1 

2 

5 


16000 

3.68(1) 

0.14(2) 

2.33(7) 

6.7 

2 

5 

5 

256000 

4.096(5) 

0.482(5) 

1.07(1) 

3.5 

1 

6 


128000 

4.082(4) 

0.482(5) 

1.07(1) 

7.4 

1 

5 


64000 

4.087(5) 

0.481(4) 

1.08(1) 

6.9 

1 

4 


32000 

4.042(7) 

0.453(6) 

1.17(2) 

7.0 

1 

4 


16000 

3.964(4) 

0.406(4) 

1.32(1) 

42 

1 

4 


256000 

4.031(9) 

0.40(1) 

1.25(3) 

1.2 

2 

7 


128000 

3.990(7) 

0.367(9) 

1.36(2) 

2.6 

2 

6 


64000 

3.950(9) 

0.34(1) 

1.46(3) 

2.1 

2 

5 


32000 

3.83(1) 

0.24(1) 

1.82(5) 

4.3 

2 

5 


16000 

3.627(8) 

0.051(9) 

2.69(5) 

30 

2 

5 

2 

128000 

4.01(1) 

0.542(9) 

1.04(3) 

0.8 

1 

5 


64000 

4.023(6) 

0.551(5) 

1.06(1) 

3.9 

1 

4 


32000 

4.01(1) 

0.54(1) 

1.09(3) 

1.6 

1 

4 


16000 

3.94(2) 

0.49(1) 

1.23(4) 

1.5 

1 

4 


128000 

3.95(2) 

0.47(2) 

1.25(5) 

0.2 

2 

5 


64000 

3.920(9) 

0.44(1) 

1.34(4) 

0.8 

2 

5 


32000 

3.84(2) 

0.35(3) 

1.58(8) 

0.4 

2 

5 


16000 

3.69(3) 

0.22(3) 

2 . 1 ( 1 ) 

1.8 

2 

5 


Table 5: The results of the fits to Eq. (|63|) for 7 ii(r;iV) for pure gravity (m = 2), 
Ising (m = 3) and three-states Potts (m = 5) model coupled to gravity. 
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m 

N t 

c 

a 

C 

dh 

x 2 

^ min 

T max 

3 

256000 

1.75(2) 

0.46(2) 

0.73(3) 

4.13(3) 

0.26 

l 

6 


128000 

1.74(2) 

0.43(2) 

0.76(2) 

4.10(2) 

0.56 

l 

5 


64000 

1.72(2) 

0.42(2) 

0.77(2) 

4.08(3) 

0.44 

l 

4 


256000 

1.68(3) 

0.30(7) 

0.87(7) 

4.01(5) 

0.26 

2 

7 


128000 

1.62(3) 

0 . 21 ( 6 ) 

0.97(6) 

3.94(4) 

0.74 

2 

6 


64000 

1.58(3) 

0.16(7) 

1.04(8) 

3.87(5) 

0.21 

2 

5 

5 

256000 

1.56(2) 

0.35(2) 

0.73(2) 

4.27(3) 

1.15 

1 

6 


128000 

1.56(1) 

0.36(2) 

0.73(2) 

4.27(2) 

0.81 

1 

5 


64000 

1.57(2) 

0.38(2) 

0.71(2) 

4.28(3) 

0.67 

1 

4 


256000 

1.43(3) 

0.07(7) 

0.97(6) 

4.06(5) 

0.46 

2 

7 


128000 

1.43(3) 

0.08(6) 

0.98(6) 

4.05(5) 

0.85 

2 

6 


64000 

1.36(3) 

-0.03(7) 

1.10(7) 

3.93(5) 

0.88 

2 

5 


Table 6 : The results of the fits to Eq. (|63| ) for n v (r]N) for Ising (m = 3) and 
three-states Potts (m = 5) model coupled to gravity. 


m 

N t 

c 

a 

C 

dh 

x 2 

7* min 

T* max 

3 

256000 

1.32(2) 

0.51(4) 

0.49(2) 

3.97(6) 

0.13 

l 

6 



1.34(1) 

0.53(3) 

0.50(1) 

4.00(4) 

0.15 

l 

9 



1.37(1) 

0.59(2) 

0.54(1) 

4.10(3) 

1.3 

l 

13 



1.38(2) 

0 . 66 ( 8 ) 

0.56(4) 

4.13(7) 

0.14 

2 

11 



1.44(2) 

0.86(7) 

0.66(4) 

4.32(6) 

1.5 

2 

14 


128000 

1.35(2) 

0.56(4) 

0.52(2) 

4.03(7) 

0.07 

1 

5 



1.38(1) 

0.61(3) 

0.55(2) 

4.14(4) 

0.47 

1 

8 



1.40(1) 

0.65(2) 

0.58(1) 

4.20(3) 

1.2 

1 

10 



1.44(3) 

0.81(8) 

0.65(4) 

4.32(8) 

0.26 

2 

9 



1.51(2) 

1.02(7) 

0.77(4) 

4.53(6) 

1.8 

2 

12 


64000 

1.35(2) 

0.54(3) 

0.51(2) 

4.05(5) 

0.15 

1 

5 



1.39(1) 

0.61(2) 

0.55(1) 

4.16(3) 

1.9 

1 

8 

5 

256000 

1.57(2) 

0.55(3) 

0.47(2) 

3.93(5) 

0.08 

1 

6 



1.59(1) 

0.59(3) 

0.50(2) 

4.00(3) 

0.54 

1 

9 



1.63(1) 

0.64(2) 

0.53(1) 

4.06(3) 

1.6 

1 

11 



1.67(3) 

0.79(8) 

0.61(4) 

4.18(7) 

0.53 

2 

10 



1.73(2) 

0.94(7) 

0.70(5) 

4.33(6) 

1.7 

2 

12 


128000 

1.55(2) 

0.53(3) 

0.46(2) 

3.86(5) 

5.1xl0 ~ 3 

1 

5 



1.57(1) 

0.57(2) 

0.48(1) 

3.93(3) 

0.67 

1 

8 



1.59(1) 

0.60(2) 

0.50(1) 

3.98(3) 

1.7 

1 

9 



1.61(3) 

0 . 68 ( 8 ) 

0.53(4) 

4.03(7) 

0.42 

2 

8 



1.65(3) 

0.77(7) 

0.58(4) 

4.12(6) 

1.1 

2 

9 


64000 

1.55(2) 

0.53(3) 

0.45(2) 

3.87(4) 

0.43 

1 

5 



1.60(1) 

0.61(2) 

0.51(1) 

4.01(3) 

2.9 

1 

7 


Table 7: The results of the fits to Eq. (| 66 | ) for the normalized spin-spin correlation 
function for Ising {m = 3) and three-states Potts (m = 5) model coupled to gravity. 
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dh 

m — 2 

m — 3 

m — 5 

Method 

4.05(15) 

4.11(10) 

4.01(9) 

ni(r; N) 

FSS 

3.92-4.01 

3.99-4.08 

3.99-4.10 

7 ii(r; N ) 

SDS 

3.85-3.98 

3.96-4.14 

4.05-4.20 

(■ l n ( r ))N 

FSS 


4.28(17) 

4.46(33) 

ntp{r\ N) 

FSS 


3.90-4.16 

4.00-4.30 

n<p(r; N) 

SDS 


3.96-4.38 

3.97-4.39 

G<p(t)n 

SDS 


Table 8: A summary of the results for dh shown in Table [ij Table |7]. FSS in the 
Method column stands for Finite Size Scaling and SDS for Small Distance Scaling. 
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Figure 1: (a) dh(a) from collapsing ni(r; N ) for Nr = 8000-64000 for pure gravity. 
Data is collapsed in groups of three lattice sizes and in the graph we indicate the 
largest of each group. We show the errors computed from x 2 only f° r the largest 
lattice in order to simplify the graph. The errors for the smaller lattices are quite 
similar. ( b ) Same as in (a) for the Ising model for N T = 8000-128000. (c) Same as 
in ( b ) for the three-states Potts model. 
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Figure 2: (a) dh(a) from collapsing n^(r;iV) for iVr 
model the same way as described in Fig. [l]. (b )Same 
Potts model. 


= 8000-128000 for the Ising 
as in (a) for the three-states 
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-2 -1 0 


log(x) 

Figure 3: ( a ) The small x behaviour of the logarithmic derivative of ni(r;N). We 
use N t = 32000-128000, dh = 4.02 and a = 0,0.54. ( b ) Same as in (a) for the 
Ising model coupled to gravity. We plot for N T = 64000-256000, a = 0, 0.51 and 
dh. = 4.08. (c) Same as in ( b ) for the three-states Potts model coupled to gravity 
where a = 0, 0.48 and dh = 4.10. 

28 







Figure 4: (a) The small x behaviour of the logarithmic derivative of n v (r;iV) for 
the Ising model coupled to gravity. We use Nt = 16000-256000 and x is obtained 
by using dh = 4.13, a = 0.45. ( b ) Same for the three-states Potts model coupled to 
gravity. We use now dh = 4.27 and a = 0.35. 
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Figure 5: (a) The rescaled according to Eq. (65) normalized spin-spin correlation 
function g 9 {r\ N) for the Ising model coupled to gravity. We use N T = 2000-256000, 
dh = 4.0, a = 0.51. ( b ) Same for the three-states Potts model coupled to gravity. 
We use a = 0.55. 
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Figure 6: (a) The small x behaviour of the logarithmic derivative of the rescaled 
normalized spin-spin correlation function N) for the Ising model coupled to 
gravity. We use N T = 16000-256000 and x is obtained by using dh = 4.0, a = 0.51. 
( b ) Same for the three-states Potts model coupled to gravity. We use now dh = 4.0 
and a = 0.55. 


31 





2 


x 


CM 

LL 


1.5 


0.5 




0 O B ^* 


2000 

4000 

8000 

16000 

32000 

64000 


•% 


0 


2 


1.5 


1 


0.5 


0 


“L 


2000 

4000 

8000 

16000 

32000 

64000 



0 


2 


x 


2 

1.5 

1 

0.5 

0 

0 1 2 


-Jr*™ 


\ 


2000 

4000 

8000 

16000 

32000 

64000 


% 




X 


Figure 7: (a) The (Z 2 (r ))at distributions for pure gravity rescaled according to 

Eq. (|60| ) using dh = 3.88, a = 0.00. ( b ) Same as in (a) for the Ising model coupled to 
gravity, x is scaled using d h = 3.99, a = 0.0. (c) Same as in ( b ) for the three-states 
Potts model coupled to gravity, x is scaled using d h = 4.07, a = 0.0. 
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Figure 8: c4(a) from collapsing (/ 2 (r ))jv for iVy = 16000-64000 for the Ising model 
coupled to gravity. 
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-2 -1 0 
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Figure 9: (o)The small x behaviour of the logarithmic derivative for pure gravity 
for (/ n (r ))n for Nt = 2000-64000. ( b ) Same as in (a) for the Ising model coupled 
to gravity and (c) for the three-states Potts model coupled to gravity. 
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Figure 10: The loop length distribution function for pure gravity for Nt = 64000. 
The dashed line is a fit to Eq. fl3ip . 
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Figure 11: (a) The loop length distribution function for the Ising model for Nt = 
64000 shown together with the pure gravity one. A wider range in y is included in 
the plot in order to show the m-dependence in the areas where finite size effects are 
important. ( b ) same as in (a) for the three-states Potts model. 
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Figure 12: The loop length distribution function for pure gravity (c = 0) and the 
c = —2 model for Nt = 64000. The data for the c = —2 model is taken from JTO . 
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Figure 13: (a)The loop length distribution function for the three-states Potts model 
coupled to gravity for Nr = 64000 and 32000 for fixed l. ( b ) The same, but for 
(approximately) fixed l/N 2 ^ dh for dh = 4.0. Similar plots are obtained for pure 
gravity and the Ising model coupled to gravity. 
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